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ELECTRON  IONIZATION  CROSS  SECTIONS  IN 
THE  DISTORTED-WAVE  APPROXIMATION 

I.  INTRODUCTION 

The  purpose  of  this  report  is  to  present  results  of  distorted-wave 
(DW)  calculations  of  electron  impact  ionization  cross  sections.  Our 
investigation  has  been  stimulated  by  the  need  for  reliable  cross  section 
data  for  a  large  number  of  ions  for  which  experimental  results  are  not 
available,  and  by  the  inaccuracy  of  most  theoretical  methods  which  often 
overestimate  the  value  of  the  ionization  cross  sections. 

The  accuracy  of  various  approximations  usually  becomes  better  as 
the  energy  of  incident  electrons  increases  and  therefore  we  have  concen¬ 
trated  our  study  on  the  low-energy  region  where  the  reliability  of  dif¬ 
ferent  methods  can  be  easily  assessed.  The  exact  quantum-mechanical 
evaluation  of  an  ionization  cross  section  is  not  possible  because  it  would 
involve  a  prohibitively  large  volume  of  numerical  calculations.  Our  goal 
was  to  develop  an  approximation  based  on  the  quantum-mechanical  approach 
to  the  collision  problem,  simple  enough  so  that  it  would  be  suitable  for 
practical  calculations,  and  such  that  it  would  give  more  reliable  results 
than  similar  methods.  We  have  investigated  several  modifications  of  the 
distorted-wave  approximation  which  appear  to  be  appropriate  for  this 
purpose. 

In  Section  II  we  present  an  outline  of  the  general  theory  of  ionization 
for  hydrogen  and  hydrogenic  ions  and  in  Section  III  we  discuss  approximate 
methods  which  will  be  used  for  comparison  with  our  own  results.  The  dis¬ 
torted-wave  approximation  without  exchange  is  described  in  Section  IV. 

We  introduce  four  versions  of  the  DW  approximation  which  are  character¬ 
ized  by  different  orthogonal  properties  of  wave  functions  and  by  different 

choices  of  asymptotic  charges  for  the  ejected  electron.  Results  of  these 
Manuscript  submitted  April  15,  1980. 
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four  modifications  are  compared  for  H,  He,  He+,  Li+,  and  Mg+.  The  beat 
approximation  is  then  used  for  the  calculation  of  ionization  cross  sec¬ 
tions  for  more  complex  atoms.  In  Section  V  we  study  the  effect  of  exchange 
on  the  ionization  of  hydrogen.  Results  are  presented  for  four  exchange 
approximations  in  which  different  choices  of  asymptotic  charges  were  used 
in  the  exchange  amplitude.  The  reliability  of  all  approximations  is  dis¬ 
cussed  in  Section  VI. 

II.  GENERAL  THEORY  FOR  HYDROCEN  AND  HYDROGENIC  IONS 

1-3 

The  total  cross  section  for  ionization  is  given  by  the  expression 

E-E. 

Q  (E)  =  i  J  o(ko,k,  k')d(ik2),  (1) 

o 

1  2 

where  E  *  j  kQ  is  the  energy  of  the  incident  electron,  E^  the  ionization 

energy,  k  and"k*  are  momenta  of  the  two  free  electrons  after  an  ionizing 

collision,  and  o(k  ,  k,  k')  is  the  single  differential  cross  section 
o 

(SDCS)  for  ionization.  All  quantities  in  (1)  and  throughout  this  paper 

are  expressed  in  atomic  units  except  the  cross  section,  which  will  be 

2  -17  2 

given  in  units  of  iraQ  *  8.7974  x  10  cm  .  Since  the  two  electrons  in 
the  final  state  cannot  be  distinguished  from  each  other, 

a(k  k,  k')  -  a (k  ,  k',  k)  (2) 

o  o 

12  1 

and  the  integral  in  (1)  is  symmetric  with  respect  to  j  k  -  y  (E-E^. 

The  integration  can  be  cut  off  at  y  (E-E^)  so  that 

r  i  2 

Q  (E)  *  J  a(kQ,  k)  d(|  k^)  (3) 

o 

with  a<k  ,  k)  -  a(k  ,  k,  k')  for  k  <  k' . 
oo  — 


2 


SDCS  for  the  ionization  of  the  ground  state  of  hydrogen  and  hydro- 

1-3 

genic  ions  can  be  written  in  terms  of  scattering  amplitudes  in  the  form 


o (kQ,  k,  k') 


k  k'  f 

ir  k  J 
o  J 


|f(£  k"')|  2  +  I g (k »  k*)|  2 


-  1 f I  | g | cos  (B  -  6') 


d  k  d  k'  . 


(4) 


f(k,  k')  and  g(k,  k’)  are  the  direct  and  the  exchange  scattering  ampli¬ 
tude,  respectively ,  k  and  k*  are  unit  vectors  in  the  direction  of  k  and 
k ' ,  and 

6  -  arg  f.  S'  *  arg  g. 

4 

It  has  been  shown  by  Peterkop  that  the  direct  and  exchange  amplitudes 
are  related  by 

f(£,  it’)  -  g(£\  k)  (5) 


so  that  (2)  is  satisfied. 


The  integral  expression  for  f(k,  k’)  is’ 


f (k,  k')  -  -  (2tt)~  a  exp 


[i  A(k,  k’)] 


/  *  r2) 


(6) 


(H  '  “total ^  *  (rl*  V  drl  dr2  • 


^^rl’  r2^  C^e  exact  unsymmetrized  solution  of  the  SchriJdinger  equation 
describing  the  scattering  of  the  incident  electron  with  the  coordinate 
r2  on  the  atom  with  the  bound  electron  which  has  the  coordinate  r^.  H  is 
the  total  Hamiltonian,  Etotaj.  is  the  energy  of  the  system  and  $  (r^,  r^) 
is  a  function  with  asymptotic  form 


*z)~  ®  (z*  r"i)  ®(z’*  ;  *2)* 

where  the  functions  'P  satisfy  equations  of  the  form 

[v2  +  k2  +  ~  +  2V(r)j  9  (z,  -k;  ^)-  0 


(7) 


(8) 


and  asymptotically  describe  incident  and  scattered  waves.  The  continuum 
functions  <p  are  normalized  per  unit  electron  density  at  infinity.  v(r)  may 
be  any  short-range  potential  and  the  effective  charges  Z  and  Z’  are  related 
by 


Z  V_ 

k  k' 


-£  ♦ 


k-k* 


(9) 


with  Zq  equal  to  the  residual  charge  of  the  atom  or  ion  after  ionization. 
Equation  (9)  follows  from  the  condition  that  the  scattering  amplitude  shall 
have  no  divergent  phase  factor.  It  has  a  simple  classical  interpretation: 
at  large  distances,  ^  v  k  t,  r2~k't,  where  t  is  the  time,  and  con¬ 
sequently  (9)  can  be  written  as 


—  +  — 
rl  r2 


Z  Z 

^  +  _o 
rl  r2 


^rrV 


The  charges  Z  and  Z’  can  therefore  be  interpreted  as  angle-dependent 
quantities  which  asymptotically  take  full  account  of  the  Coulomb  potentials. 
The  phase  A(k,  k')  in  (6)  is  given  by 


A  (k 


2[f  -(*)♦  £  -  (f) 


(10) 


with  4  X  *»  E  -  E^,  provided  that  (9)  is  satisfied. 
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III.  APPROXIMATE  METHODS 

Approximate  methods  for  the  evaluation  of  the  scattering  amplitude 
(6)  must  be  used  because  the  exact  solution  ^Cr^,  r^)  of  the  Schrddinger 
equation  is  not  known.  Moreover,  it  is  difficult  to  satisfy  condition 
(9)  for  effective  charges  Z  and  Z’  as  they  depend  on  the  angles  between 
k  and  k'.  If  (9)  is  not  satisfied,  then  the  relative  phases  of  the  direct 
and  exchange  amplitude  remain  undetermined  and  have  to  be  prescribed  in¬ 
dependently  of  any  other  approximation. 

Approximation  methods  fall  in  two  categories:  non-exchange  approxi¬ 
mations  in  which  the  exchange  amplitude  gain  (4)  is  neglected,  and  approxi¬ 
mations  where  exchange  is  taken  fully  into  account .  Non-exchange  approxi¬ 
mations  are  computationally  much  simpler  and  may  be  successfully  used  for 
high  energies  of  incident  electrons,  when  exchange  is  relatively  unimportant. 
The  effect  of  exchange  may  sometimes  be  reduced  by  partial  cancellation  of 
the  two  terms  in  (4)  containing  g,  which  explains  why  the  non-exchange 
approximations  often  give  satisfactory  results  for  complex  atoms. 

Our  discussion  of  various  approximate  methods  will  be  restricted  to 
a  few  basic  approximations  which  will  be  used  later  for  comparison  with 
our  final  results . 

3.1  Non-exchange  Approximations. 

These  approximations  assume  that  g  *  0  in  equation  (4).  The  value 
of  the  phase  45^,"k')  does  not  affect  the  result. 

The  simplest  and  most  widely  used  method  is  the  Bom  approximation 
for  neutral  atoms.  For  the  ionization  of  hydrogen  from  the  Is  state  it 
adopts  the  conventions 
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Z  -  Z  ,  Zf  =  Z  -  1, 
o  o  ’ 

Y  (r1>  r2)  -  rx)  exp  (ikQ  -r2), 

(11) 

$  (^1»  •  cp(Zq,  exP  C-ik^'*  r^)  * 

V(r)  -  0. 

In  the  Coulomb-Born  (CB)  approximation  for  hydrogenic  ions  with  nuclear 

charge  Z  »  exp  (ilT  •  rl)  in  (11)  is  replaced  by  cp  (Z  -1,  k*  ;  rt.)  and 
00/  o  o  z 

exp  C-ikT'  •  r^)  by  cp(Z' ,  -  k' ;  r^). 

In  the  non-exchange  approximations,  the  two  electrons  in  the  final 

state  are  considered  to  be  distinguishable  and  therefore  (2)  is  not  valid. 

Consequently  the  expression  (1)  for  the  total  cross  section  is  not 

identical  to  (3).  Using  (l)for  the  calculation  of  the  total  cross  section, 

one  has  to  omit  the  factor  *s  because  it  is  associated  with  the  indistin- 

guishability  of  electrons.  Then  (1)  together  with  (11)  gives  the  Born  (a) 

3 

approximation  (according  to  the  notation  introduced  by  Rudge  ) .  On  the 
other  hand,  equation  (3)  leads  to  the  Born  (b)  approximation. 

In  the  Born (a)  method,  the  ejected  electron  (with  the  coordinate 
r^)  is  associated  with  the  effective  charge  Zq  for  all  values  of  the 
momentum  k.  This  has  an  unsatisfactory  consequence  for  k  >k',  because 
the  effective  charge  acting  on  the  faster  electron  is  larger  than  that 
acting  on  the  slower  one.  A  comparison  of  the  Bom  (a)  and  Bom(b) 
approximations  for  hydrogen  is  shown  on  Fig.  1.  The  Bom(b)  method 
usually  gives  substantially  better  results  than  the  Bom  (a)  approximation. 
Both  the  Bom  and  the  CB  methods  can  be  easily  generalized  to  raany-electron 
atoms  and  ions. 
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1  14-20  42 

Many  results  from  these  approximations  are  now  available  ’  ’  * 

44  45 

’  .  Various  improvements  have  been  attempted  involving  better  approxi¬ 

mations  to  Y  and  different  choices  of  Z  and  Z'.  The  distortion  of  wave 
functions  describing  the  free  electrons  in  both  the  initial  and  final 
states  can  be  taken  into  account  by  properly  defining  V(r)  in  Eq.  (8). 
Approximations  with  V(r)  ^  0  are  usually  referred  to  as  distorted-wave 
(DW)  approximations. 

3.2  Exchange  Approximations 

The  exchange  amplitude  g(k,  kT' )  can  be  calculated  from  (5),  and 
if  Eq.  (9)  is  not  satisfied,  the  value  of  cos  (8-8')  has  to  be  specified. 

By  setting  cos  (8-8')  *  0  and  using  (11),  we  obtain  an  approximation  which 
is  identical  to  the  Bom (a)  method.  With  the  choice  cos  (8~8')  =  1  we 
get  the  Bom-exchange  approximation  due  to  Peterkop^.  This  choice 
corresponds  to  maximum  interference  and  therefore  gives  the  smallest  cross 
section  of  any  approximation  for  the  scattering  amplitude.  The  Bom- 
exchange  approximation  represents  a  substantial  improvement  in  the  calcu¬ 
lated  cross  sections.  It  gives  excellent  agreement  with  experiment  for  the 
ionization  of  the  ground  state  of  He  ,  but  for  H  it  still  overestimates  the 
cross  sections  by  about  30%  at  low  energies  (Fig.  3). 

Another  way  to  obtain  the  exchange  amplitude  is  to  use  expression 
(6)  for  the  direct  amplitude  and  to  replace  <$(r^,  r^)  by  $(r^,  r^).  This 
procedure,  together  with  the  choice  cos  (8-8')  ■  1  and  Eq.  (3)  leads  to 
the  Born-Oppenheiraer  (BO)  approximation.  An  attractive  feature  of  this 
method  is  the  possibility  of  deriving  both  the  direct  and  exchange 
scattering  amplitude  from  the  expression  analogous  to  (6)  using 
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antisymmetric  function  4*  and  taking  into  account  the  spin  of  the  electrons  . 
In  the  BO  approximation,  (5)  is  satisfied  if  for  k  >  k'  in  the  direct 
amplitude  we  define 


$(r^»  r”2)  =  exp  ("ik*"  *^L)  cp(Zo,  -k'  ;  r*).  (12) 

Then  Eq.  (1)  is  identical  to  (3).  As  a  consequence,  the  smaller  effective 

charge  is  always  associated  with  the  faster  of  the  two  electrons,  but  the 

effective  charges  for  individual  electrons  are  discontinuous  at  k  =  k'. 

Another  disadvantage  of  the  BO  method  is  the  lack  of  orthogonality  of  $ 

defined  by  (12)  and  Y  given  by  (11).  For  the  ionization  of  hydrogen,  the 

14 

BO  approximation  grossly  overestimates  cross  sections  at  lower  energies. 

The  possibility  of  various  choices  of  phases  and  effective  charges 

associated  with  the  direct  and  exchange  amplitude  leads  to  a  large  variety 
3  6  7  21  23 

of  approximations  ’  ’  ’  .  However,  for  the  ionization  of  hydrogen, 

no  method  has  given  results  substantially  better  than  the  Bom-exchange 
approximation . 

IV.  Ionization  Cross  Sections  in  the  Non-Exchange  DW  Approximation 

In  this  section  we  outline  the  approximation  used  in  the  present 
paper  for  the  calculation  of  ionization  cross  sections  without  electron 
exchange.  Four  different  approximations  (A1  -  A4)  will  be  investigated 
and  results  will  be  presented  for  H,  He,  He+,  Li+,  and  Mg+.  Approximation 
A4  which  appears  to  be  most  reliable  will  then  be  used  for  the  calculation 
of  ionizaflon  cross  sections  of  more  complex  atomic  systems. 
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4.1  General  Description 

In  our  approximation,  distorted  waves  are  used  to  describe  the 
incident,  scattered,  and  ejected  electrons.  For  the  ionization  of 
hydrogen  and  hydrogenic  ions  we  set 


r2)  -  iKls;  r^)  «(ZQ-1,  k^;  r^)  , 

(13) 

*(rlf 

r* )  -  cp(Z ,  -lT;  r^)  tp(Z',  -k ' ;  r2). 

(14) 

The  short-range  potential  V(r)  corresponding  to  cp  describing  the  incident 
electrons  in  (13)  represents  an  additional  effect  of  nuclear  charge  and 
the  screening  of  the  bound  electron.  We  put 


V(r)  ■  vls(r)»  with 


(15) 


Vr> 


rrf  l2 

rm  r 

J  P(ls ;  x)J  dx  + 

1  X_1  [p(ls;  x)_ 

dx.  (16) 


Leaving  the  choice  of  effective  charges  Z  and  Z’  for  the  time  being  open, 
we  define  V(r)  corresponding  to  the  functions  ®  in  (14)  by 


V (r)  =>  0(Z,  Zq)  (Zq  -  Z)  Vls(r) 


(17) 


with  a  similar  expression  for  Z'  instead  of  Z.  S(Z,  Z  )  is  the  step 
function  given  by 

e(Z,  Z  )  =  0  for  Z  >  Z  ,  0(Z,  Z  )  =  1  for  Z  <  Z  .  (18) 

O  0  0  —  o 

Definition  (17)  formally  guarantees  a  gradual  change  from  the  full 

screening  by  the  bound  electron  to  no  screening  if  Zq  *  Z  or  Z'. 
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Ionization  cross  sections  are  calculated  from  Eq.  (3)  and  (4)  with 

g  *  0.  We  now  generalize  this  approximation  to  complex  atoms  and  ions. 

Let  us  consider  a  collision  process  in  which  an  electron  is  ejected  from 

the  n£^q  shell.  To  obtain  the  expression  for  the  single  differential 

cross  section  o,  we  have  to  sum  over  all  possible  final  states  of  the 

ejected  electron  and  over  all  final  terms  n£^q  ^  S'L'  of  the  remaining 

ion.  Using  partial  wave  expansions  for  the  functions  c o  in  (13)  and  (14) 

2  26 

and  assuming  LS  coupling,  we  obtain  ’ 


cr(n£  q;  k  ,  k)  ■=  16q  k"2  (2A+1)"1  (2£+l)  (2£'+l)  (21'  +1) 

1  °  °  1 


C  :■ ft V  [«*.  ■;  ■■)] 


(19) 


where  2.^,  I,  2.j  and  2,'  are  angular  momentum  quantum  numbers  of  the  bound, 
incident,  ejected  (slower),  and  scattered  (faster)  electron,  respectively. 
Parameter  A  corresponds  to  the  multipole  order  in  the  expansion  for  the 
interaction  potential.  The  direct  radial  integral  is  given  by 


RX  (V;  ll  V)  "  J  P(n)ll;  rl)  FCZo-i.  k0»  l’  r2} 


«(X,  0)  U(rx) 


rX  ■* 

+  6  (A ,0)  U'(r2)+ 

r  _ 


(20) 


F(Z,  k,  rx)  F(Z|  k|  l';  r2)  d^  dr2 


with 


U(rx) 


(Z  -  Z  ) 
o 


+  0(Z, 


V 


(21) 
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p  .ipmifiippii 


U'(r2)  -  (Z1  -  Zo)  [  rf  +  e(Z',Z0)  (tj)]  .  (22) 

P(n£^;  r)  is  the  radial  function  of  the  bound  electron  in  the  n£^  orbit 
and  the  continuum  functions  F(a,  b,  c;  r)  are  solutions  of 

[4-  — f— -  +  ^7  +  2V(r)  +  b2]  F(a,  b,  c;  r)  -0  (23) 

L  dr  r  c 


normalized  so  that  at  infinity 


F(a,  b,  c;  r)  ~  b  **  sin  £br  +  ab  ^  An  (2cr)  +  nj  • 


(24) 


(r)  is  given  by  (16)  if  P(ls;  x)  is  replaced  by  P(ni^;  x) .  For 
complex  atoms  we  set 

V (r)  =  V  (r)  +  8  (A,  Z  )  (Z  -  A)  V  (r) ,  (25) 

c  o  o 

where  Z  is  the  residual  charge  of  the  ion  after  ionization,  A  is  equal  to 
o 

Z  or  Z'  for  the  continuum  functions  of  the  final  state  and  to  Z  -  1  for 

o 

the  initial  state,  8 (A,  Z  )  -is  defined  by  (18),  and 

o 


Vc(r) 


Nr 


{r_1/  [vx)]  dx  +  /  x_1  [Pi(x)] 

i®l  o  r 


dx 


(26) 


The  summation  in  (26)  is  carried  out  over  all  bound  electrons  in  the  final 
state  described  by  radial  functions  P^.  In  this  procedure  we  take  into 
account  only  the  spherical  component  of  the  interaction  potential  in 
calculating  the  distortion  of  continuum  wave  functions. 

Unlike  in  the  Born(b)  and  CB(b)  methods,  the  wave  functions  of  the 
initial  and  final  state  of  the  total  system  with  distorted  continuum 
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functions  representing  free  electrons  are  generally  not  orthogonal.  In 
our  DW  approximation,  the  orthogonality  is  achieved  by  solving  Eq.  (23) 
for  F(Z,  k,  i  r^)  with  the  right-hand  side  replaced  by  PCnJl^;  r^)  if 
£|  *  ^l"  Parameter  is  adjusted  so  that  J  P(n£^;  r)  F(Z,  k,  2,^;  r)  dr 
*  0.  Then  expression  (21)  for  U(r^)  will  contain  an  additional  term 

P(n4i;  rL),  (27) 


! 


but  radial  integrals  (20)  will  not  depend  on  U'^)  due  to  the  orthogonality 
of  P(n£^;  r)  and  F(Z,  k,  r). 

If  the  initial  or  final  state  of  the  system  is  represented  by  an 
antisymmetric  function,  the  procedure  described  above  is  not  sufficient 
to  assure  orthogonality.  This  may  be  achieved  by  orthogonalizing  also  the 
continuum  function  F(Zj  k',  l*  ;  r)  of  the  scattered  (faster)  electron  to 
P(n£^;  r) .  Although  such  an  additional  orthogonalization  is  not  required 
by  necessity  in  non-exchange  approximations,  it  may  change  the  computed 
cross  sections  and  we  will  investigate  its  effects  in  the  following 
sections.  We  will  also  Introduce  two  different  choices  of  effective 
charge  Z,  thus  having  four  approximations,  A1  -  A4. 
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4.2  Approximation  Al. 

In  order  to  keep  our  calculations  as  simple  as  possible,  we  will 
assume  that  charges  Z  and  Z'  are  independent  of  the  angles  between  k  and 
k ' .  This  automatically  implies  that  Eq.  (9)  will  not  be  satisfied  and 
we  must  use  other  considerations  to  determine  their  values. 

First  of  all,  the  violation  of  condition  (9)  may  not  be  as  serious 


as  it  appears,  because  in  the  calculation  of  scattering  amplitudes 


according  to  (6),  the  most  important  parts  of  the  cp  functions  are  those 
close  to  the  atom.  In  this  region,  the  form  of  cp  is  predominantly  deter¬ 
mined  by  all  possible  interactions  of  atomic  and  free  electrons,  of  which 
only  some  were  taken  into  account  by  the  form  of  V(r)  in  (8).  By  a  proper 
adjustment  of  Z  and  Z*  one  may  include  effects  of  other  neglected  inter¬ 
actions  and  achieve  a  better  representation  of  the  functions  cp  in  this 
important  region,  even  if  (9)  is  not  satisfied. 

Secondly,  Eq.  (9)  is  not  sufficient  by  itself  to  determine  uniquely 
both  Z  and  Z',  and  if  Z'  is  arbitrarily  chosen,  then  (9)  may  lead  to  a 
physically  unsatisfactory  value  of  Z.  For  example,  if  we  choose  Z'  «=  0  for 
neutral  atoms  and  if  both  electrons  move  in  the  same  direction,  then  from 

(9)  Z  *  1  + - -  <  1  and  therefore  the  corresponding  cd  will  des- 

k'  k'~k  -2  -2 
cribe  the  slower  electron  as  being  subject  to  the  force  -Z  r^  , 

-2  -2 

while  actually  it  is  subject  to  the  force  -r^  -  <  -r.^2 

In  the  Born  and  the  CB  approximations  it  is  assumed  that  the  full 
residual  charge  Zq  is  acting  on  the  slower  electron  in  the  final  state. 

Eq.  (9)  then  takes  the  form 

Z’  -  Zq  -  k’  \k  -  k'  I"1  (28) 

and  the  choice  Z'  ■  Zq  -  1  is  compatible  with  the  spherical  average  of 
(28)  over  all  directions  of  k. 

In  the  A1  approximation  we  adopt  the  same  choice  of  Z  and  Z'  and  set 

Z  -  Z  ,  Z'  -  Z  -  1.  (29) 

o  o 

Due  to  the  value  of  Z' ,  the  continuum  functions  F(Z',  k',  for  the 

scattered  electron  are  orthogonal  to  F(Z  -1,  k  ,  l;  r«)  and  therefore  the 

o  o  ^ 

radial  integrals  (20)  are  independent  of  U(r^)  given  by  (21)  and  (27). 


jf 


4.3  Approximation  A2 

The  choice  of  Z  and  Z'  is  the  same  as  in  Al,  but  we  use  continuum 
functions  F(Z',  k' ,  l';  r)  orthogonalized  to  P(n2^;  r)  if  2'  *  2^.  This 
procedure  introduces  an  additional  term 

-  e2  P(nilS  r2)  (30) 

into  the  expression  (22)  for  U'Gtj)*  but  that  does  not  affect  radial 
integrals  (20)  as  they  do  not  depend  on  U’^).  On  the  other  hand,  the 
radial  integral  with  2  *  2’  *  2^  m  i.  '  will  now  contain  a  non-vanishing 
contribution  corresponding  to  U(r^)  due  to  the  lack  or  orthogonality  of 
F(Zo~l,  kQ,  2;  r and  F(Z',  k’,  2,';  r~)  -  This  contribution  is  presumably 
small  and  will  be  ignored  in  the  present  approximation. 

4.4  Approximation  A3 

In  considering  the  choice  of  Z  and  Z’  for  the  direct  scattering 
amplitude,  we  start  from  the  experimental  observation  that,  for  incident 
energies  not  far  above  the  threshold,  the  double  differential  cross 
section  (DDCS)  for  ionization  is  larger  for  scattering  angles  t?  <  90° 

‘chan  for  &  >  90°.  This  applies  to  both  the  faster  and  the  slower  elec¬ 
trons.  For  the  faster  electrons,  the  DDCS  has  a  strong  peak  in  the  forward 
direction,  while  the  DDCS  for  slower  electrons  shows  a  much  less  pronounced 
maximum  around  t?  ■  60°.  This  can  be  interpreted  as  a  manifestation  of  the 
fact  that  both  electrons  move  predominantly  in  the  forward  direction, 
especially  when  they  are  close  together,  before  their  mutual  repulsion 
increases  the  angle  between  their  velocities.  Therefore,  instead  of  taking 
a  spherical  average  of  |k  -  k'|  in  (9),  we  determine  Z  and  Z'  by  assuming 
that  k  is  parallel  to  k The  effective  charges  may  be  interpreted  in 
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several  ways  and  we  obtain  various  conditions  for  Z  and  Z’: 

(a)  In  accordance  with  (9)  and  our  assumption  about 
k  and  k  ,  we  get 


Z_ 

r. 


r2"rl 


(31) 


using  r^  ~  kt,  r^  ~  k't, 

(b)  If  we  require  the  actual  force  acting  on  the  slower 

2 

electron  to  be  -  — j,  we  obtain 

rl 


(r2"rl^ 


(32) 


(c)  The  same  condition  written  for  the  faster  electron  is 

2 


Z’ 

l 

r2 


(r2-rx) 


(33) 


(d)  Classically  both  electrons  move  in  such  a  way  that  the  force 

Z  Z 

acting  on  the  center  of  gravity  is  equal  to - j - • 

rl  r2 

We  may  require  the  same  force  to  act  on  the  center  of 
gravity  if  the  electrons  are  moving  in  the  field  of  effective 
charges  Z  and  Z'  and  we  obtain 


Z’ 


(34) 


All  four  conditions  are  simultaneously  satisfied  if  Z  and  Z'  are  chosen 
according  to  (32)  and  (33).  For  r^/r2  +0  we  obtain 
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r 


Z  -  Zo  +  (Vr2)‘ 


z’~  zo  -  1  -  cri/r2r 


(.35) 


However,  our  basic  assumption  about  the  parallelism  of  k  and  k'  cannot 

be  valid  in  the  limit  of  r^/ r2  -*-1,  and  conditions  (31)  -  (33)  are  not 

applicable.  On  the  other  hand,  Eq.  (34)  does  not  contain  the  term 

(r2~r^)  ^  and  is  therefore  more  appropriate  for  the  determination  of 

Z  and  Z'  in  the  limit  of  r^/r2  =  1*  Generally,  if  lT  is  not  parallel  to 

k',  the  radial  component  of  the  force  acting  on  the  center  of  gravity  of 

Z  Z 

the  two  electrons  is  equal  to--^-  cos  ct - cos  g  t  where  a  and  8  are 

rl  r2 

angles  between  the  direction  to  the  center  or  gravity  and  the  vectors 
k  and  k  ,  respectively.  If  the  effective  charges  produce  the  same  force, 
then. 


cos  a  - 


Z' 


cos  6 


Z  Z 

O  0  0 

— j  cos  a - j  cos  8  . 


(36) 


For  small  values  of  r-^/ r2 ,  (36)  is  equivalent  to  (34)  due  to  the  assumption 
of  parallelism  between  k  and  k'.  When  r^/r2  -*-1,  a  approaches  8  regardless 
of  the  angle  between  k  and  k’  and  in  the  limit  again  (36)  becomes  equivalent 
to  (34).  Therefore  Eq.  (34)  is  approximately  satisfied  for  all  values  of 
r^/r2.  Thus  we  adopt  the  following  expression  for  Z: 


Z  -  Z 

o 


(37) 


If  the  scattering  amplitude  is  calculated  from  (6)  using  functions 
(13)  and  (14)  ,  it  is  much  less  sensitive  to  the  value  of  Z’  (associated 
with  the  faster  electron)  than  to  Z.  In  order  to  simplify  the  calcu¬ 
lations,  we  therefore  adopt  a  constant  value  of  Z'  and  define 
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Z’  -  Zo  -  1  (38) 

in  agreement  with  (35)  in  the  limit  of  r^r ^  “  0.  This  choice  is  consistent 
with  the  Bora  and  CB  methods  and  with  our  approximations  A1  and  A 2. 

Putting  r^/r ^  ~  k/k',  Eq.  (37)  together  with  (38)  gives  finally 

Z  -  Z  +  (k/k’)2.  (39) 

o 

Eq.  (39)  agrees  with  (35)  which  is  valid  for  small  values  of  k/k'. 

According  to  this  choice,  the  attractive  force  acting  on  the  slower  elec¬ 
tron  is  !arger  than  the  force  caused  by  the  residual  charge  of  the  ion. 

As  a  result,  it  is  more  difficult  for  the  ejected  electron  to  escape  and 
this  leads  to  the  lowering  of  the  ionization  cross  section.  An  unsatis¬ 
factory  property  of  the  choice  (38)  and  (39)  is  that  in  the  limit  of 
k  *  k'  the  effective  charges  of  the  two  electrons  differ  by  2.  Moreover, 
the  assumption  that  both  electrons  move  predominantly  in  the  forward 
direction  is  of  questionable  validity  near  the  threshold  of  ionization. 

Similarly  to  the  A1  approximation,  the  continuum  functions  for  the 
incident  and  scattered  electrons  are  orthogonal  and  all  terms  in  radial 
integrals  associated  with  U(r^)  and  U'^)  vanish. 

4 . 5  Approximation  A4 

The  choice  of  Z  and  Z'  is  the  same  as  in  A3  and  is  given  by  (38) 
and  (39).  In  A4,  continuum  functions  F(Z',  k',  l';  r)  with  V  *  are 
orthogonalized  to  P(n£^;  r)  as  in  A2  and  resulting  additional  terms  in 
the  radial  integral  associated  with  U(r^)  are  ignored. 

17 

. -  ^  ..  ...t.  V.""' 'j !« 


Basic  characteristics  of  approximations  A1  -  A4  together  with  the 
Born  and  CB  methods  are  summarized  in  Table  1. 


Table  1 


Approximation 

Z 

Bom  and  CB 

Z 

o 

A1 

Z 

o 

A2 

Z 

0 

A3 

Z  +(k/k')2 

0 

A4 

Z  +(k/k')2 
o 

Function  of  the 

Scattered  (slower)  I 

electron 

Orthogonal  to  P(ni,^; 

;  r) 

Orthogonal  to  P(nX,^; 

;  r) 

4.6  Results  and  Discussion 

In  the  numerical  evaluation  of  (19)  we  have  included  the  first  five 

angular  momenta  corresponding  to  the  ejected  (slower)  electron.  That 

was  found  sufficient  in  most  cases  in  the  low-energy  region  with  E/E^  <_  4. 

The  summation  over  l  and  V  has  been  carried  out  until  convergence  has 

been  reached.  The  number  of  angular  momenta  taken  into  account  was 

2 

usually  less  than  20.  For  each  value  of  E  3  's  kQ  we  have  computed 

a(k  ,  k)  for  4  values  of  k.  Hartree-Fock  wave  functions  of  bound  electrons 
o 

29 

calculated  by  Clement!  and  Roetti  have  been  used  for  non-hydrogenic 
atomic  systems. 

The  differences  between  various  approximations  diminish  as  the 


incident  electron  energy  increases.  At  the  same  time  the  number  of 
angular  momenta  l  and  £'  which  should  be  taken  into  account  becomes 


prohibitively  large.  Therefore  we  have  made  calculations  in  the  low- 
energy  region,  where  the  reliability  of  different  methods  can  be  easily 
assessed. 

Results  from  approximations  A1  -  A4  for  H,  He+,  He,  Li+,  and  Mg+ 
are  shown  in  Table  2.  In  Figs.l  -  5  they  are  compared  with  experimental 
data  and  with  other  approximations.  Two  general  conclusions  can  be 
drawn  from  Table  2:  Firstly,  the  use  of  continuum  functions  for  the 

scattered  (slower)  electron  orthogonalized  to  the  bound  wave  function 
(approximations  A2  and  A4)  always  reduces  the  value  of  the  cross  section, 
and  secondly,  the  choice  of  effective  charge  Z  given  by  (39)  in  A3  and 
A4  has  the  same  effect.  The  differences  are  most  pronounced  in  the  low- 
energy  region  as  expected,  and  the  effect  of  orthogonalization  is  much 
smaller  for  ions  than  for  neutral  atoms.  Approximation  A4  yields  the 
lowest  cross  sections  and  appears  to  be  the  best  in  comparison  with 
experiment. 

For  hydrogen  (Fig.  1),  A1  and  A2  give  results  similar  to  the  Bom(b) 
approximation  except  near  the  threshold.  Approximations  A3  and  A4  lead 
to  a  substantial  improvement  of  results,  especially  in  the  energy  region 
where  the  cross  section  reaches  its  maximum  value.  Very  close  to  the 
threshold,  A4  appears  to  underestimate  the  cross  sections  for  H  and  He. 
However,  in  the  case  of  hydrogen,  the  value  0.076  for  E/E^^  =  1.25  represents 
the  lower  limit,  because  omitted  contributions  from  >  4  in  (19)  may 
not  be  negligible.  Measurements  for  hydrogen  represented  by  curves  3 
and  4  in  Fig.  1  were  normalized  to  the  Born  approximation  at  high  energy. 

If  these  measurements  were  normalized  to  absolute  data  of  Rothe  et  al.11 
shown  by  curve  5,  the  maximum  value  of  the  cross  section  would  be  reduced 


by  about  0.1  f  a  .  As  a  result,  even  the  A4  approximation  would  still 
slightly  overestimate  the  cross  section  in  the  maximum  region,  but  the 
agreement  in  the  low-energy  region  would  be  improved. 

Approximations  A3  and  A4  agree  very  well  with  experimental  data  for 

He+  at  all  energies  (Fig.  2).  For  the  Helium-like  systems  (Figs.  3  and 

4)  there  is  a  discrepancy  in  the  region  of  maximum  cross  section.  A4 

overestimates  the  value  of  Q  by  about  20%  and  the  difference  is  probably 

due  to  the  omission  of  exchange,  which  is  known  to  reduce  substantially 

32 

excitation  cross  sections  in  helium-like  ions 

For  Mg  ,  A1  and  A2  give  very  large  cross  sections  compared  to 
experiment  and  also  to  the  CB(b)  approximation  (Fig.  5).  A  possible 
explanation  for  this  difference  may  be  the  focusing  effect  of  the  short- 
range  potential  V(r)  on  incident  electrons  in  the  DW  approximation.  This 
effect  may  be  compensated  in  A3  and  A4  by  a  larger  attractive  force 
acting  on  the  ejected  electron  so  that  a  good  agreement  with  experimental 
data  is  achieved. 

In  Table  3  we  present  ionization  cross  sections  calculated  in  approxi¬ 
mation  A4  for  several  atoms  and  ions  of  increasing  complexity.  The  results 
for  two  representative  energies  are  compared  with  experiment  and  with 
other  theoretical  methods.  For  atoms  and  ions  with  ground  configurations 
consisting  of  several  terms  we  have  used  ionization  energies  which 
correspond  to  the  lowest  term.  It  was  also  assumed  that  all  final  terms 
have  the  same  energy  in  agreement  with  the  derivation  of  Eq.  (19). 

Effects  of  autoionization  have  been  neglected  and  the  cross  sections  for 
atoms  with  the  ns^np^  configurations  have  been  evaluated  from  both  the 
ns^  and  np^  subshells. 


The  approximation  A4  gives  smaller  ionization  cross  sections  than 

the  Coulomb-Born  method  for  almost  all  ions  shown  in  Table  3  and  generally 

agrees  better  with  experiment.  The  theoretical  cross  sections  computed 

in  A4  for  C+,  N+,  and  0+  at  E/E^  =  4  are  in  a  good  agreement  with  measured 

values.  High  experimental  cross  sections  in  the  low-energy  region  with 

E/E^  =*  2  for  these  ions  are  most  likely  to  be  attributed  to  the  effect  of 

autcionization  of  metastable  levels  which  are  populated  by  excitation 

of  the  2s  electron.  This  effect  has  been  discussed  by  Moores.^  In  the 

case  of  ions  with  closed  or  almost  closed  up  shells  (Ne+,  Na+,  K+) ,  A4 

overestimates  the  cross  sections  by  about  25%,  but  the  results  for  K+ 

43 

agree  with  measurements  of  Peart  and  Dolder.  The  worst  agreement  is 
+3 

obtained  for  C  ,  where  the  values  shown  in  Table  3  represent  a  lower 

limit  due  to  the  omission  of  terms  with  >  4.  It  is  interesting  to  note 

that  no  other  approximation  was  able  to  reproduce  the  experimental  results 

+4 

for  this  ion,  although  the  agreement  for  the  next  isoelectronic  ion,  N  , 

40 

was  found  satisfactory.  Theoretical  cross  sections  for  neutral  oxygen 
are  in  satisfactory  agreement  with  experiment  regarding  the  uncertainty 
of  measured  data. 

We  conclude  that  the  approximation  A4  is  the  best  of  all  four  non¬ 
exchange  approximations  A1  -  A4 ,  and  it  generally  yields  better  results 
than  the  Bom(b)  or  Coulomb-Bom (b)  method  except,  perhaps,  in  the  low- 
energy  region  close  to  the  ionization  threshold,  where  it  may  under¬ 
estimate  the  cross  sections. 
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V.  EFFECT  OF  EXCHANGE  IN  THE  DW  APPROXIMATION 


5.1  General  Considerations 

In  the  following  sections  we  will  investigate  the  effect  of  exchange 
scattering  amplitude  g(k,  k')  on  the  value  of  ionization  cross  sections 
and  our  discussion  will  be  restricted  to  the  ground  state  ionization  of 
hydrogen . 

The  total  cross  section  for  ionization  is  calculated  from  (3)  and 
(4)  and  the  exchange  scattering  amplitude  is  obtained  from  the  direct 
amplitude  by  employing  relation  (5).  According  to  Eq.  (3),  k  <.  k'. 

We  adopt  the  convention  that  effective  charge  Z  is  always  associated  with 
momentum  k,  i.e.  with  the  slower  electron  in  the  final  state. 

Expanding  the  continuum  functions  into  partial  waves  we  obtain  for 
the  ionization  of  H(ls) 


(k  ,  k)  =  16  k-2  V*  (2£+l)  (2£'+l)  (2£ 

o  O  faW 


( l  V  4'  \ 

;«>L  oV 


44 '  5  | 


X  j  Rd  (2X£;  4|4')J  +  ^A'+l)'1  (24^+1)  R®,  (4^5  4£« ' )J 

* 

-  (2£'+l)_1  (24’+l)  |rJ,  (414;  4j4')|  |R®,  (4j4;  4[4')|  cos  (n' )  | 


(40) 


Direct  radial  integrals  Rd  are  given  by  (20)  and  exchange  integrals 
R  are  defined  by 

A 
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R®  (AjA;  AjA')  -  j  P(ls ;  rj  F(Zq-1,  kQ,  A;  r2)  [«(X,0)  U(r2) 

(41) 

r  X 

+  6  (A  ,0)  U'(r1)  +  — ^-JF(Z,k,A^;  r2)  FUJkjA'jrj.)  d^  dr2- 

The  value  of  y-y '  is  associated  with  relative  phases  of  the  direct  and 
exchange  part  of  the  element  of  the  scattering  matrix  which  connects 
channels  characterized  by  kQ,  l  and  k,  k',  Our  choice  of  Z  and 

Z'  does  not  satisfy  condition  (9)  and  therefore  there  is  no  point  in 
attempting  to  determine  y-y'  from  Eq.  (10).  In  our  approximation  we  set 

cos  (y-y')  =  1 

in  analogy  to  the  choice  of  relative  phases  in  the  Born-exchange  approxi¬ 
mation^.  Expression  (40)  is  valid  also  for  ionization  of  hydrogen  and 
hydrogenic  ions  from  any  ns  state  if  P(ls;  r)  in  (41)  is  replaced  by 
P(ns;  r) ,  and  for  ionization  of  atoms  with  one  ns  electron  outside  closed 
shells  if  V(r)  is  properly  defined  according  to  (25)  and  (26). 

The  results  of  non-exchange  DW  approximations  presented  in  Section 
4.6  indicate  that  the  approximation  A4  gives  the  best  overall  agreement 
with  experiment.  We  therefore  adopt  the  same  choice  of  Z  and  Z'  as  in 
A4  for  the  calculation  of  direct  integrals  in  (40),  i.e. 

Z  -  Z  +  (k/k')2,  Z*  =■  Z  -  1.  (42) 

o  o 

Moreover,  in  agreement  with  A4  we  orthogonalize  both  continuum  functions 
F(Z,k,£^;  r)  and  F(Z',k',£';  r)  with  £|  -  0  and  £'  *  0  to  the  bound 
function  P(ls;  r)  as  described  in  section  4.1.  The  orthogonalization  of 
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continuum  functions  appearing  in  the  exchange  integral  is  performed  as 
well.  Due  to  the  orthogonal  property  of  P(ls;  r)  and  F(Z' ,k' ;  r) , 
the  contribution  to  the  integral  (41)  associated  with  U^)  vanishes. 

According  to  Eq.  (5)  the  exchange  amplitude  g(k,  k')  with  k  <  k' 
is  equal  to  the  direct  scattering  amplitude  for  the  process  in  which  the 
ejected  electron  has  a  higher  velocity  than  the  scattered  one.  The 
direct  amplitude  may  be  interpreted  on  the  basis  of  a  classical  picture 
of  the  scattering  process  and  one  can  see  that  the  effective  charges  Z 
and  Z'  in  the  case  when  the  ejected  electron  is  faster  need  not  neces¬ 
sarily  be  associated  with  the  charges  Z  and  Z'  when  it  is  slower.  There¬ 
fore,  even  if  Z  and  Z'  were  specified  by  (42)  for  direct  radial  integrals, 
we  have  to  determine  their  values  separately  for  exchange  integrals.  We 
will  investigate  four  different  choices  described  in  the  following 
sections . 


5 . 2  Approximation  A5 

In  the  expression  for  exchange  radial  integrals  we  put 

Z-Z  -1,  Z’-Z  +  (k/k')2  .  (43) 

o  o 

This  choice  is  analogous  to  the  effective  charges  adopted  in  the  Born- 
exchange  approximation  in  the  sense  that  the  exchange  amplitude  can  be 
obtained  from  the  direct  amplitude  by  interchanging  lT  and  k'  in  (14). 

If  the  exchange  amplitude  is  interpreted  as  a  direct  amplitude  according 
to  (5),  then  effective  charges  for  both  electrons  are  continuous  at 
k  *  k'.  This  is  shown  graphically  on  Fig.  6,  where  effective  charges 
are  plotted  against  k/k'  (with  k  <.  k')  in  such  a  way  that  the  diagram 
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covers  the  whole  interval  of  possible  velocities  of  the  ejected  electron. 
With  respect  to  the  choice  of  Z,  the  integral 

.  Jf(zo-i,  V  r>  F<z>  k>  r>  dr 

0 

in  either  zero  or  has  a  small  value  so  that  the  contribution  to 
corresponding  to  U'(r^)  may  be  neglected.  An  unsatisfactory  property  of 
the  choice  (43)  is  associated  with  the  fact  that  the  charge  for  the  slower 
electron  in  the  exchange  amplitude  is  smaller  than  the  charge  acting  on 
the  faster  electron,  a  situation  which  is  similar  to  the  Born-exchange 
approximation . 

5 . 3  Approximation  A6 

In  this  approximation,  the  exchange  amplitude  is  obtained  from 
the  direct  amplitude  by  interchanging  r^  and  r ^  in  (14)  as  in  the 
Born-Oppenheimer  approximation.  This  procedure  leads  to  the  choice 

Z  =  Z  +  (k/k')2,  Z'  =  Z  -  1  (44) 

o  o 

for  effective  charges  in  (41).  The  higher  charge  is  always  associated 
with  the  slower  of  the  two  electrons,  but  the  charges  are  discontinuous 
at  k  •  k'  (Fig.  6). 

The  functions  F(Z  -1,  k  ,  J£;  r)  and  F(Z,  k,  l';  r)  with  l  *  £' 
oo  l  i 

are  not  orthogonal  because  Z  ?  Zq-  1  and  consequently  the  contribution 
to  the  exchange  integral  corresponding  to  U'(r^)  does  not  vanish.  Never¬ 
theless,  we  will  ignore  this  contribution  in  the  present  approximation. 
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5.4  Approximation  A7 


In  analogy  to  our  assumption  about  the  parallelism  of  k  and  k' 

in  the  direct  amplitude,  we  will  assume  that,  if  the  ejected  electron  has 

->■ 

a  higher  velocity,  the  vectors  k  and  k'  are  antiparallel.  In  the  classi¬ 
cal  model  it  means  that  the  ejected  electron  moves  in  the  forward  direction, 
while  the  scattered  electron  moves  backwards.  This  assumption  takes 
explicitly  into  account  the  possibility  of  the  backscattering  which  is 
well  documented  by  experimental  evidence.  Using  the  same  arguments  as 
in  Sec.  4.4,  we  define 

Z  »  Z  -  (k/k')2,  Z'  =  Z  -  1  (45) 

o  o 

in  the  exchange  amplitude. 

In  this  approximation,  the  charge  for  electron  with  coordinate 
is  continuous  at  k  *  k'  (Fig.  6)  and  the  slower  of  both  electrons  is 
always  associated  with  a  higher  effective  charge.  An  unsatisfactory 
consequence  of  (45)  is  that  it  leads  to  full  screening  of  both  electrons 
at  k  *  k'.  "As  in  A6,  we  will  ignore  the  non-vanishing  contributions  to 
associated  with  U'(r^)  in  (41). 


5.5  Appro xima t ion  A8 


With  the  choice 


Z  -  Z'  -  Z 


adopted  in  this  approximation  for  exchange  amplitude,  U'^)  vanishes 
according  to  (22)  and  effective  charges  are  discontinuous  at  k  *  k' 

(Fig.  6).  This  choice  is  similar  to  the  one  us-^d  by  Geltman  and  Hidalgo*1 
in  the  Coulomb-projected  Born  approximation. 
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Various  characteristics  of  exchange  approximations  A5  -  A8  are 
summarized  in  Table  4. 

5.6  Results  and  Discussion 

Results  from  approximations  A5  -  A8  for  the  ground  state  ionization 
of  hydrogen  are  shown  in  Table  5  and  Fig.  7,  where  they  are  compared  with 
experiment ,  the  Born-exchange  method  and  the  non-exchange  approximation 
A4. 

The  symbol  >  indicates  that  the  number  of  angular  momenta  in¬ 
cluded  in  the  calculation  was  not  sufficient  and  therefore  the  value  shown 
in  the  table  represents  a  lower  limit.  Approximation  A5  gives  the  lowest 
cross  sections  and  the  best  agreement  with  experimental  data  if  one  takes 
into  account  the  possibility  that  a  different  normalization  according  to 
absolute  measurements  of  Rothe  et  al.1^  would  lower  curve  2  in  Fig.  7. 
However,  A5  underestimates  cross  sections  near  the  threshold,  where  better 
results  are  obtained  from  A7.  Approximation  A6  gives  very  unsatisfactory 
results  for- all  energies  in  spite  of  its  theoretical  attractiveness.  Com¬ 
pared  to  the  non-exchange  approximation  A4,  only  A5  gives  lower  cross 
sections  at  all  energies,  but  the  difference  is  not  too  pronounced  and  it 
is  smaller  than  the  difference  between  the  Born(b)  and  Born-exchange 
approximations . 

The  differences  between  various  approximations  are  more  conspicuous 

if  one  compares  the  partial  contributions  to  the  single  differential  cross 

section  a(kQ,  k).  Table  6  shows  a  typical  example  for  one  incident 

2 

energy  E  »  2Z^.  The  dominant  contribution  at  low  k  corresponds  to 

2 

«  1  in  all  approximations.  With  increasing  k  the  contributions  from 


higher  values  of  become  more  important,  notably  in  A6,  A7  and  A8. 

2 

The  importance  of  exchange  increases  with  increasing  k  and  the  dif¬ 
ferences  between  various  approximations  are  most  significant  in  the  limit 
k  *  k'.  They  also  increase  with  decreasing  energy  E.  This  is  shown  on 
Figs.  8  -  10  which  represent  the  single  differential  cross  section 
a(kQ,  k)  for  three  energies  of  incident  electrons.  The  effects  of 
exchange  are  least  pronounced  in  A5  and  one  can  see  why  the  total  cross 
sections  in  this  approximation  are  little  different  from  the  non-exchange 

approximation  A4:  the  decrease  of  a  caused  by  exchange  at  low  values 

2  2 
of  k  is  partially  compensated  by  an  increase  for  k  approaching  E  -  E^. 

VI.  SUMMARY  AND  CONCLUSION 

In  the  present  Investigation  we  have  studied  eight  different  modifi¬ 
cations  of  the  distorted-wave  approximation  for  the  calculation  of  elec¬ 
tron  ionization  cross  sections.  From  the  non-exchange  approximations 
A1  -  A4,  A4  appears  to  be  the  most  satisfactory  and  generally  gives  a 
better  agreement  with  experiment  than  the  Born  or  Coulomb-Bora  approxi¬ 
mations,  without  excessive  requirements  on  the  computing  time.  It  is 
therefore  suitable  for  a  production  of  large  number  of  data  needed  in  the 
analysis  and  interpretation  of  hot  plasmas  in  laboratory  devices  and  in 
astronomical  objects. 

Exchange  approximations  A3  -  A8  have  been  tested  on  hydrogen,  and 
A5  has  been  found  superior  to  A6  -  A8.  Ionization  cross  sections  obtained 
in  A5  are  slightly  lower  than  those  from  A4,  but  the  difference  does  not 
appear  to  be  too  significant  with  respect  to  the  uncertainty  of  experi¬ 
mental  data  used  for  comparison.  The  small  effect  of  exchange  in  the  A5 
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approximation  may  be  just  a  particular  property  of  hydrogen,  but 
good  agreement  of  the  non-exchange  approximation  A4  with  measured  cross 
sections  in  many  other  cases  suggests  more  general  validity  of  this 
property. 

For  incident  energies  corresponding  to  E/E^  >  4,  the  importance 
of  exchange  rapidly  decreases  and  in  this  region  A4  is  expected  to 
give  reliable  cross  sections  for  all  atoms  and  ions. 
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Table  2 


Ionization  Cross  Sections  (in7raQ)  in  the  Non-Exchange  DW  Approximation 


Atom 

Approxi¬ 

mation 

E/Ei 

1.25 

2 

4 

10 

H(ls) 

A1 

0.401 

0.815 

1.012 

0.629 

A2 

0.112 

0.672 

0.994 

0.629 

A3 

0.291 

0.594 

0.771 

A4 

>0.076 

0.499 

0.759 

0.535 

He+(ls) 

A1 

0.0230 

0.0577 

0.0644 

0.0394 

A2 

0.0207 

0.0563 

0.0644 

0.0394 

A3 

0.0194 

0.0487 

0.0563 

A4 

0.0173 

0.0474 

0.0563 

He  (Is2) 

A1 

0.115 

0.379 

0.556 

A2 

0.045 

0.326 

0.548 

A3 

0.120 

0.326 

0.475 

A4 

0.040 

0.283 

0.469 

+  2 

Li  (Is  ; 

A1 

0.0191 

0.0541 

0.0661 

A2 

0.0172 

0.0509 

0.0659 

A3 

0.0174 

0.0484 

0.0597 

A4 

0.0156 

0.0471 

0.0586 

hg+(3s) 

A1 

mmm 

0.788 

m 

A 2 

0.782 

B 

A3 

0.279 

0.562 

0.566 

A4 

0.269 

0.556 

0.566 

— 

Table  3 


Comparison  of  the  Non-Exchange  DW  Method  (Approximation  A4)  with  Other 
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TABLE  5 


Ionization  Cross  Sections  (in  TraQ)  for  Hydrogen 


in  the  DW-Exchange  Approximation 


Approximation 

1.25 

E/E. 

2.0 

4.0 

A5 

0.072 

0.453 

0.669 

A6 

>  0.245 

>  0.898 

>  0.820 

A7 

0.143 

0.602 

0.724 

A8 

0.340 

0.575 

0.673 

TABLE  6 

2 

Partial  Contributions  to  k  a(k  ,k)  from  angular  momenta  l!  of  the 
_ _ o _ o _ _ _ i _ 

Slower  Electron  for  E/E^  ■  2.  Ionization  of  H(ls) 


Approxi¬ 

mation 

2 

k 

l ' 

1 

0 

1 

2 

3 

4 

A4 

0 

1.021 

7.783 

1.228 

0.0636 

0.0014 

0.1667 

0.352 

3.217 

1.185 

0.193 

0.1176 

0.3333 

0.0485 

0.979 

0.780 

0.260 

0.0459 

0.5 

0.0415 

0.112 

0.346 

0.234 

0.0738 

A5 

0 

1.021 

7.783 

ESOi 

0.0636 

0.0014 

0.1667 

0.270 

2.559 

Ell  -its: 

0.172 

0.0158 

0.3333 

0.0934 

0.885 

0.625 

0.213 

0.0379 

0.5 

0.0661 

0.475 

0.385 

0.212 

0.0661 

A6 

0 

1.088 

6.672 

1.384 

0.113 

0.0044 

0.1667 

0.471 

3.058 

1.202 

1.933 

0.0172 

0.3333 

0.380 

1.546 

2.451 

1.643 

0.671 

0.5 

0.432 

0.679 

1.843 

2.600 

2.030 

A7 

0 

1.088 

6.672 

1 

Egg 

0.0044 

0.1667 

0.525 

3.300 

ESS 

0.0355 

0.3333 

0.413 

1.929 

0.956 

0.304 

0.0706 

0.5 

0.303 

0.882 

0.440 

0.221 

0.0741 

A8 

0 

1.060 

6.161 

1.071 

0.0622 

0.0018 

0.1667 

0.492 

2.711 

1.148 

0.210 

0.0251 

0.3333 

0.327 

1.513 

1.125 

0.390 

0.0886 

0.5 

0.437 

1.635 

1.358 

0.668 

0.229 
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measurement  of  Peart  and  Dolder33 ;  other  symbols  as  for  Fig.  1. 


LOS(E/Ej ) 

Fig.  5  —  Cross  sections  for  ionization  of  Mg+(3s).  1  —  Coulomb-Born  (b)l-®,  ejection  of  3s 
electron;  2  —  Coulomb-Bom  (b)^®;ejection  of  2p  electron;  3-measurement  of  Martin  et  al. 31- 
other  symbols  as  for  Figure  1. 


DIRECT  AMPLITUDE  EXCHANGE  AMPLITUDE 


f 

T 

S 

Fig.  6  —  Effective  charges  acting  on  electron  1  (full  line) 
and  2  (broken  line)  in  approximations  A5-A8.  s  —  slower 
electron  (charge  Z),  f  —  faster  electron  (charge  Z') 


mu 


If 


sections  for  the  ground  state  ionization  of  hydrogen.  1  —  Bom  exchanged 
2  —  measurement  of  Boksenberg9.10;  3  —  A4 
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Fig.  8  —  Single  differential  cross  section  a(k0,  k)  for  hydrogen. 
E/E;  *  1.25.  Full  line  —  A4;  broken  lines  —  exchange  approxi¬ 
mations  A5-  A8.  Omitted  terms  in  A4  and  A6  would  increase 
a  for  k2  close  to  E  —  Ej. 


8.  E/Ei  *  2.  Omitted  terms  in  A6 
ise  a  for  k2  close  to  E  —  Es. 


